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The stability of an evaporating thin Hquid film on a solid substrate is investigated within 
lubrication theory. The heat flux due to evaporation induces thermal gradients; the gener- 
ated Marangoni stresses are accounted for. Assuming the gas phase at rest, the dynamics 
of the vapour reduces to diffusion. The boundary condition at the interface couples trans- 
fer from the liquid to its vapour and diffusion flux. The evolution of the film is governed 
by a lubrication equation coupled with the Laplace problem associated with quasi-static 
diffusion. The linear stability of a flat film is studied in this general framework. The 
subsequent analysis is restricted to diffusion-limited evaporation for which the gas phase 
is saturated in vapour in the vicinity of the interface. The stability depends then only on 
two control parameters, the capillary and Marangoni numbers. The Marangoni effect is 
destabilising whereas capillarity and evaporation are stabilising pr ocesses. The results o f 
the linear stability analysis are compared with the experiments of lPoulard et ol\ l)200,'^ 
performed in a different geometry. In order to study the resulting patterns, an amplitude 
equation is obtained through a systematic multiple-scale expansion. The evaporation rate 
is needed and is computed perturbatively by solving the Laplace problem for the diffusion 
of vapour. The bifurcation from the flat state is found to be a supercritical transition. 
Moreover, it appears that the non-local nature of the diffusion problem unusually affects 
the amplitude equation. 



1. Introduction 

Since the pioneering studies of lThomsorJ ill 85,'^ . lMa,ra,ngonil l)l 8fi,'j) a,nd lRena,rdl l|l 9nft . 

much attention has been devoted to the now called Marangoni instabilities. Thomson and 
Marangoni first proposed surface tension gradients as a cause for convection in liquids. 
The Marangoni effect consists in the variation of surface tension with temperature or 
liquid composition and drives this class of instabilities. The hexagonal patterns observed 
by^cnard (1900) in thin layers heated from b elow prompted a number of studies (for 
reviews, see lDaviall987tlSchatz fc Neitze]|l200l|) . Recent research in this field has focused 
on the correct description of the gas above the fluid layer and of the deformability of 
the interface ( VanHook et al. 19 97; Golovin et al, .1997), or th e effect of local heating 
l)Miladinova et a^Jl2003 iKalliadasis et a/Jl2003t lYeo a]jl26ol . 

Marangoni effect can als o be driven by e vaporation. Many experimental situations 
of interest are reviewed bv iBerg. Boudartfc Acrivos (1966). On the one hand, evap- 
oration generates thermal gradients as the phase transformation requires latent heat. 
For sp reading droplets of slightly volatile liquids, lR,edon. Bro chard-Wvart fc Rondeled 
l|l992l) reported festoon instabilities near the contact line while .Kavelipour et ah (,2002^") 
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measured height fluctuations over the whole drop. lHegseth. R.a shidnia' IZ^haHigge') ob- 
served vigorous interior flow in evaporating droplets of a volatile liquid. On the other 
hand, in the case of mixtures, a difference of the evaporation rate between compo- 
nents changes the relative conc entrations at the interface and so generates surface ten- 
sion gradients. iFburnier fc Cazabat (1992 ) : Vuilleumier. Ego. Neltner fc Cazabat (199^; 
iFanton &: CazabatI lll998l):lHosoi fc ijushTllBOOll) studied tears of wine and the associated 
convection rolls. ^ NguvenfcSteb^ l)2002|) showed instabilities induced by surfactants. This 
compos ition mechanism may enhance evaporation and is therefore usefu l in drying tech- 
niques l|Marra fc Huethorslill99it IO'Brienlll993 iMatar fc Crasteill200l|) . 

A comprehensive theoretical stability analysis of evaporating/condensing films was 
done by B.ur£lb.ach. Bankoff & Davis (,1958.) . T hey included vapour recoil and thermo- 
capillarity. but as in sub sequent studies (see lOron et all Il997t iMargerit et 
iMerkt fc BestehornI l2003l) . the evaporation and condensation are governed by the de- 
parture from thermodynamic equilibrium at the in terface. Within this framewor k evapo- 
ration is intrinsically destabilising as can be seen in lProsperetti fc Plessetl ( Il984l) who did 
not consider Marangoni stresses. In the case of very thin films ('Elbaum & Lipson 19941) 

jscj 



microscopic forces may be destabilising as sh own bv.Samid-M erzel. Lipson fc T annhauscij 
l)l998|) : lLvushnin. Golovin fc PismenI l)20o3) . It is worth noticing here that all models de- 
veloped in these papers are one-sided: they do not account for the gas phase dynamics 
except throu gh the boundary condition at the interface. 

In contrast . iDeegan et al\ l)l997lll99^ showed that evaporation of pinned water droplets 
is limited by diffusion of vapour in air, thus they explained the origin of coffee stains, 
[^^achile et ah (,2002) explained their experiments on freely receding evaporating droplets 
within the same framework. They also observed the drops of certain fluids to loose their 
axisymmetry and the contact line to become wavy (|Poulard et a/.l |2003tl. These unex- 
plained instabilities are among the motivations of the present study. In particular, the 
established one-sided model would always predict an instability for evaporat ion. 

Our aim here is to generalize the one-sided study of'Burelb ach et al\ l(l988l) for evapo- 
rating thin films by taking into account the dynamics of the vapour. In section|21 we build 
a model which includes both thermodynamically determined transfer of the molecules 
across the interface and diffusion of vapour in the gas phase. This generalises the two 
class of models presented above. We describe the liquid film within lubrication theory 
taking into account surface tension gradients and loss of mass. In section 3, we perform 
the full linear stability analysis of this system. Then we restrict to the diffusion limited 
regime which is relevant for the experiments of ^Poulard et al. (2003J. We find the ampli- 
tude equation describing patterns above the in stability onset. Event ually, in section 4, 
we compare our results with the experiments of lPoulard et al\ l)2003|) . 



2. The model 

We consider the dynamics of a two-dimensional bi-layered liquid-gas system over a 
solid substrate (figure 0). The gas phase is a mixture of an inert gas and of the vapour of 
the liquid which is volatile. We assume that the gas phase is not saturated by the vapour 
so that the liquid evaporates. The typical corresponding experimental situation is that 
of a water layer evaporating in air. The latent heat needed for the phase transformation 
drives a heat flux to the interface in the liquid. The induced temperature variations may 
generate surface tension gradients. 

The model derived below is built on the lubrication approximation for the liquid layer, 
and accounts for surface tension variations and loss of mass through evaporation. The 
gas phase is at rest; its dynamics is reduced to the diffusion of the vapour in the mixture. 
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Figure 1 . Geometry of the physical system. 



For a complete description of the system, the boundary conditions at the interface are 
needed and, in particular, the evaporation rate must be specified. Although this analysis 
is focused on evaporating films, it applies to the case of condensation as well. 

The physical situation of interest is that of volatile liquids such as water and alkanes 
evaporating in air. A compilation of the physical parameters is given in table ^ The 
governing equations will be obtained within some approximations, which are relevant to 
these liquids and will be justified at the end of this section. 

We restrict the study to a two dimensional system which coordinates are x along the 
substrate and z in the normal direction (figure^. The state of the system is determined 
by the height h{x, t) of the interface and the density z, t) of vapour in the gas phase 
(both are functions of space and time). In the following, we write the equations describing 
the flow and heat diffusion in the liquid phase, diffusion in the gas phase and then the 
boundary conditions at the interface. 

2.1. The liquid film 

We consider a very thin film so that we work within the long-wavelength approximation 
(lubrication theory) where the typical height H is much smaller than the typical hori- 
zontal scale and we neglect gravit y. Our starting point for the ev olution of the thin film 
is a lubrication equation following lOron. Davis fc 



dh d ( h^ 
dt dx \ 3/i 



■ d^h ^ d /J' 



2/i 9x j" pt ^ ^ 



where fi is the viscosity of the fluid and 7 is the surface tension. The last term in the 
brackets comes from the shear stress due to surface tension gradients. The right-hand 
side corresponds to the volume loss through evaporation, is the liquid density, pv the 
vapour density and J the mass flux across the interface. Let Jq be the characteristic 
evaporation rate. The ratio = Jo/ Pi is the relevant velocity scale in the system. 
Vapour thrust comes from the back reaction of molecules leaving the liquid into the gas 
phase; it is a quadratic effect in the evaporation rate. Other forces such as Van der Waals 
are not included. 

2.2. Surface tension gradient 

To close equation (I2.1|l we first need to compute the surface tension gradient. We use the 
standard linear equation of state 

7(r)-7o-|^l(r- (2.2) 
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where 70 = 7(T'subs)- This approximation is accurate on a large temperature range for 
most common hquids (not too close from a phase transition). We use the absolute value 
of = Tsubs) as surface tension decreases with temperature for most liquids. The 

temperature field T(x, z) satisfies the standard convection-diffusion equation in the liquid 
which reduces to = in the long-wavelength approximation. Neglecting density, 
viscosity, thermal diffusivity and kinetic ener gy of the gas the energy balance at the 
interface gives the heat flux ijOron et a/.lll997j) 

dT 

K-^iz — h(x)) — —CyJ 
az 

Cv being the vaporisation latent heat per unit mass, k the thermal conductivity of the 
liquid. We obtain 

T{z = h{x)) = ^^2d^h{x) + T.ubs (2.3) 

K 

where the substrate is assumed isothermal. Equations p. 212. 3(1 result in 

d7 



7 = 70 + 



dT 



^hJ. (2.4) 



To close the system (|2.1I2.4|) one has to compute the evaporation rate J. 

It will prove useful in the following to find an upper bound for the amplitude of 
the temperature variation on the interface. With the help of ((2.3(1 . we obtain: AT = 
CvJoAh/n, Ah being the height fluctuations. The height fluctuations Ah are smaller 
than H. Thus, the reduced temperature is bounded by 6* = ^%~t^i ^ being 

the characteristic thickness of the film. From table 2, we see that 6 is very small for the 
liquids we are interested in. 



2.3. The vapour 

The gas phase is at rest, so that there is only diffusion of the vapour. We consider the 
limit of quasi-static diffusion where the characteristic diffusion time is much smaller than 
the characteristic evaporation time H^/D <C Hpi/Jq, D being the diffusion coefficient 
of the vapour in the gas phase. In terms of the Peclet number Pe = fov§, the latter 
condition is Pe ^ 1. Hence, the vapour concentration p{x, z, t) (local number of particles 
per unit volume in the gas phase) is a solution of Laplace's equation: 

VV = (2.5) 

^ ~ '§x^ ^ '§z^ being the 2D-Laplacian. 

The gas phase is not saturated by the vapour. This condition is enforced by a constant 
diffusion rate at infinity 

-T. 7T' (2-6) 

oz D 

Expe rimentally, either the gas is pumped at t he top of the containe r i Mancini fc Mazal 
l2nf)3|) or the temperature of a top plate is fixed ( VanHook et a^llOflTi) . In both situations, 
the gas density is imposed at a certain height above the film, which induces a density 
gradient. Here, we impose the value of this gradient assuming that the height of the gas 
phase is much larger than other lengths in the system. 

To solve ((2.511 a boundary condition at the interface is needed. It is obtained in the 
next subsection. 
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2.4. Evaporation rate 

The vapour and the hquid are coupled through the evaporation rate. The kinetic theory 
leads to a linear constitutive relation between the mass a nd the departure from equilib - 
rium at the interface, known as the Hertz-Knudsen relation fProsperetti fc Plessetl l|l984) . 



^moi - '^{pl^iT,^,) - p|,„,)n (2.7) 

where M is the molecular weight, p®'^ is the density of the gas at the liquid/gas co- 
existence, p|int = p(z = h{x)) is the gas density at the interface, fee is the Boltzmann 
constant, a is the accommodation coefficient (close to unity) and n is the outward normal 

to the interface. We note uth — ''2JIT ' which is a typical kinetic velocity. 

In the gas phase, the vapour mass flux, related to the departure from uniform vapour 
density is given by: 

J = -DVp. (2.8) 
Due to the continuity of the normal evaporative flux at the interface, we have 

-D{n ■ V)p|i„t = x'th(p°^(Tint) - Plint). (2.9) 

Writing a linear temperature dependant equation of state and using H2.3(l . one obtains 
for the equilibrium density at the interface: 

p^<i(T)=p^^(T.,b.)-^^W. 

dJ K 

Thus, the boundary condition at the interface H2.9I) may be rewritten as: 

[l + Vt^^^^h^ (n • V)pU = I'th(p^'»(Teubs) - Plint). (2.10) 

This is a mixed boundary condition in the sense that it relates the value of the density 
and its normal gradient at the interface. It corresponds to the conservation of the mass 
i.e the equality of the evaporation rate and the mass flux in the gas. It includes both 
diffusion and transfer across the interface which were used separately in the litterature 
(see subsections 2.6, 2.7 and 4.1). 

2.5. Governing equations. Non dimensional parameters 

For the rescaling of the equations, we choose the typical thickness H of the liquid layer, 
the characteristic evaporation time H/vcv = Hpi/Jo, the mass flux far from the substrate 
Jo and the JqH/D as units, respectively, of length, time, evaporation rate and density 
of vapour in the gas phase. We make the substitutions h — > Hh, x — > Hx^ z Hz, 
J PiV^vJ, p p°HT^uhs) + p{JoH/D) in equations i'lAllAl and i}M'2.6V2.U)\i . The 
lubrication equation becomes 

where we have introduced Ca = Ma = Hv„vPe — — \^\ and A = """""^ respec- 

tively the capillary, the Marangoni and the vapour thrust numbers. The Laplace problem 
reads 
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Water 


Nonane 


Octane 


Heptane 


Hexane 


pi (kg/m^) 


1000 


717 


699 


682 


656 


£v (J/kg) 
7 (kg/s^) 


2.4 10'' 


3.18 10^ 


3.82 10^ 


3.21 10^ 


3.22 lO'^ 


7.20 10"^ 


22.38 10"=* 


21.77 10"^ 


20.31 10"^ 


18.42 10"^ 


M (kg/m/s) 


8.9 10"^ 


6.65 10"* 


5.08 10"* 


3.87 10"* 


3.0 10"* 


^ (kg/mVK) 
ll^l (kg/sVK) 


1.4 10"^ 


1.7 10"^ 


5.3 10"^ 


1.1 10"^ 


2.7 10"^ 


1.5 10"'' 


9.4 10"^ 


9.53 10"^ 


9.80 10"^ 


1.02 10~* 


Uev (m/s) 


0.48 10"" 


0.2 10"" 


0.6 10"" 


2.2 10"" 


6.8 10"" 


Uth (m/s) 


148 


52 


55 


58.8 


63.3 


K (kg.m/sVK) 


0.60 


1.31 10"^ 


1.28 10"^ 


1.23 10"^ 


1.20 10"^ 


fair (m^/s) 


1.9 10"^ 


3 10"" 


3 10"" 


3 10"" 


3 10"" 



Table 1. Values of physical parameters for different fluids at P = 1 atm and T = 25° C 
l|Lidelll995ll : liquid density pt, latent heat £v, surface tension 7, liquid viscosity p., temperature 
derivative of gas density at equilibrium temperature derivative of surface tension 

evaporation velocity iicv, thermal (kinetic) velocity Uth, thermal conductivity of the liquid k 
and diffusion coefficie nt of the v apour in air Dair- Wcv is estimated from the experiments on 
evaporating drops Poulard et al\ OOOSl using the formula Wcv — jo/R where jo is the evaporation 
parameter and R is the drop radius; the diffusion coefficient in air Dair is roughly estimated 
using kinetic theory; using Clapeyron relation, we have = ^yS^ ~ ^ where Ai; is the 
volume change associated which the v aporizationand p""^ is evaluated using the vapour pressure 
at saturation under 1 atm l|Lidelll99,'J) . 



with the two boundary conditions 

(l + x/i)(n-V)p|i„t = Pekp|int and lim |? = -1. (2.13) 

(jz 

where Pek — ^^^jy- is a kinetic Peclet number and x — Hvn^ ^ is called the thermal 
expansion number. 

The evaporation rate is given by 

J^-{n-V)pU (2.14) 
which is the non dimensional version of (|2.8|l . 

2.6. The one-sided limit 

The evolution equation given by Burel bach et "dl lll98?j) : IOron et al\ lll997l) can be re- 
covered as a particular limit of our model. The case of no diffusion corresponds to the 
limit of small Pek. Introducing a new scaling for the density p = pPe\./x, we find that p 
satisfies the Laplace equation (|2.12|) with the boundary conditions 

XJ(1 + X^) = -Plint and Jim ^ (2.15) 

In the limit Pek 0, the boundary condition at infinity yields a uniform density distri- 
bution p = Cst. The lubrication equation reads 

dh J_d_(l^d_\d^ A f XP V1 ^.h"^ 9 ( Xhp \ \ _ XP 

-dt ~^ c^diXYdi ['di^ [YT^J \ " 2di[YTxh) j - TT^ ^ ' 

where J has been eliminated usi ng (|2.15(l . Equation H2.16|) is the same as the equation ob- 
tained bv lBiirelbach et all 1 1988 1 up to a choice in scalings, when omitting the disjoining 
pressure term. 
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Number 
Ca (xlO**) 
Ma (xW^) 
Pe (xlO^) 

X 

X 

e (xio®) 

A (xlO") 



H 



Definition 

1 

CyVcv Pi I d7 
Kj I dT I 

V S- 



D{dp«q/dT)£v 



Pv70 

Table 2. 



Water Nonane Octane Heptane Hexane 



Signification 

viscous stresses 
capillary stresses 
Marangoni stresses 
capillary stresses 
evaporation time 
ditlusipn time 
kinetic time 
ditlusion time 
density fluctuations 
diffusion limited v.s. 
one-sided regime 
temperature fluctuations 

vapour thrust 
capillarity 

Typical values for the non-dimensional parameters. 



0.6 


0.6 


1.4 


4.2 


11 


0.7 


0.3 


1.1 


3.8 


13 


4.9 


2.0 


6.1 


22 


69 


1.5 


3.5 


3.7 


3.9 


4.2 


0.16 


0.05 


0.20 


0.38 


1.03 


9.4 


70 


18.5 


10.3 


4.1 


1.3 


0.4 


1.1 


4.0 


12 


2.8 


0.6 


1.9 


9.0 


31 



2.7. The pure diffusion limit 

We now consider the opposite limit Pe^ — > oo which we refer to as diffusion limited 
regime. The boundary conditions H2.13(l for the Laplace problem H2.12|l become 

p|int-0 and lim ^ = -I. (2.17) 

It appears that the gas is saturated in vapour immediately above the interface (in 
dimensional quantities, p\int — P^^iTsuhs)) s^nd that evaporation is limite d by diffusion. 
This bounda r y condition wa s used i n the study of evapor ating droplets bv lDeeean et all 
l|l997Lll999|) : IPoulard etai] l)2003|) : ICachile et all l)2002|) . The Laplace problem has an 
electrostatic equivalent: the one of finding the electric potential (p) with an imposed 
electric field (Jo) at infinity and a fixed constant potential on a deformed plane. The 
sharp edge effect implies a larger evaporation rate at crests which tends to restore the 
flat state so evaporation is a stabilising mechanism in the diffusion limited regime. 



2.8. Discussion 

The evaluation of the relevant dimensionless parameters for water and different alkanes 
(tableEJ shows that the limits Pe — > 0, — > and A — > are reasonable. The smallness of 
the Peclet number Pe ensures that the time needed to build up the concentration profile 
above the film is much smaller than the characteristic time for the motion of the inter- 
face, so that stationary diffusion is a good approximation. The smallness of the reduced 
temperature 9 allows a linear approximation for the gas density at the interface. The 
relevant physics is therefore contained in the values of the capillary, the Marangoni, the 
kinetic Peclet and the thermal expansion numbers (Ca, Ma, Pek, x)- Other mechanisms, 
such as molecular interactions, are neglected. 

The closed system to be studied consists in the lubrication equation H2.11|l coupled 
to a Laplace problem H2.12I2.13|I . This unusual coupling comes from evaporation which 
relates the film mass loss to the gradient of the vapour concentration. This induces non- 
locality in the lubrication equation as the mass loss is a function of the whole shape 
of the interface. To simplify notations, we drop from now on the tildes for the rescaled 
variables. 
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3. Stability of the flat interface 

Equations H2.11I2.12|) have as solution for the film thickness h{x,t) — 1 — t and gas 
density p{x, z,t) = —z + C — {l — x/Pek)t with C — {Pck — 1 — x}/P^k- As this base 
state is non stationary, linearisation of the equations gives a non autonomous partial 
differential equation, so that standard linear stability (modal) analysis should not apply. 
For simplicity, we assume from now on that the base state is h{x, t) = 1 and p{x, z, t) = 
—z + C, which amounts to adding a volume source Vov in the right-hand side of equation 
1)2.1(1 , as this source compensates exactly for the loss of mass at infinity. This also amounts 
to a quasi-steady approximation for which evaporation is sufficiently slow so that the 
thickness of the layer remains approximately constant during the growth of unstable 
modes. 



3.1. Full linear stability analysis 

We study the stability of the flat state by seeking solutions of equations H2.11I2.14() in 
the form: 

h=l + Sh, p = C - z + 6p, J =1 + SJ. 

After linearisation, those equations admit Fourier-mode solutions of wavenumber k and 
growth rate Q: 6h = A e^* Qikx_^ ^ stands for the complex conjugate of the preceding 
term); we now compute the corresponding Sp and SJ. To begin with, as a harmonic 
function, 6p has to be of the form B e**"'^ e"'*^'^ -I- c.c (B is a complex valued function of k), 
where the absolute value must be taken in order to ensure vanishing Sp at z +oo. The 
boundary condition H2.13() at the interface gives at linear order Sp\z=i — ^i^^f^^^^iJ^p^^ Sh. 
Hence, we can compute B{k) 

= „,/^:^, ^fee-l^-K-^) + c.c. 
\k\{x + l) + Pek 

From H2.14|l we get the perturbed evaporation rate 5J = — ^<5p|z=i- Plugging in H2.11|l 
and dropping non-linear contributions gives the dispersion relation: 

r2(fc) is growth rate of the wavenumber k; if Q{k) > 0, then the perturbation grows and 
the corresponding mode is unstable. Evaporation can either stabilise or destabilise large 
wavelengths depending on the values of and x whereas capillarity stabilises short 
wavelengths; Marangoni effect drives the instability. 

3.2. Transfer rate limited regime: linear stability analysis 

In the limit Pet <C min{x, fc, kx}, corresponding to the one-sided model, the dispersion 
relation 1)3.1(1 becomes 

m-~^k^ + ^k' + T^- (3.2) 

3Ca 2Ca 1 + X 

meani ng that the film is always unstable. This limit was studied in detail by bv lBurelbach et al\ 
l|l988j) . In fact, it is a singular limit (the base state has a different form), and spurious 
effects arise at fc = : equation 1(3.2(1 predicts exponential growth of a constant change 
in the film thickness. A ccounting for the time dependance of the base state as done by 
iBurelbach et all l(l988l) corrects this artifact (it yields in particular ri(0) = 0) but does 
not change the unstable behaviour. 
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Figure 2. DifFusion- limited regime: growth rate Q of the mode of wave number k (dimensionless 
quantities) for three typical value of the Marangoni number Ma. We have chosen here a capillary 
number Ca = 10~^. 



3.3. Diffusion limited regime: linear stability analysis 

We now consider the opposite limit Pe^ ^ max{x, k} which is reasonable given the 
experimental parameters. In this regime the dispersion relation ()3.1|l becomes 

The last term shows again that evaporation stabilises long wavelengths in this regime 
(figme|2Jl. The absolute value of the wavenumber k comes from quasi-static diffusion. Such 
a non-a nalyticitv in th e dispersion relation is well-known in the context of diffusive growth 
(see e.g. lLangeilll98Cl|) . The relevant control parameters are the capillary and Marangoni 
numbers. Suppose that Ca is fixed; when the Marangoni number Ma is small, there is 
no unstable mode. The first unstable wave number fee (marginally stable mode) appears 
when the Marangoni number reaches a critical value Mac such that f7(fc — kc) — and 
^17(fc = kc) = 0. Solving this system, we obtain 

Mac = 181/3 Ca2/3 

Mac (3.4) 



k'^ 



for the critical parameters at the threshold. 

The relevance of these results to the experiments is examined in the last section. 



3.4. Diffusion limited regime: weakly non-linear analysis 

The preceding linear stability analysis shows that the system may become unstable, it 
gives the most unstable wavelength and the instability threshold. However, predicting 
the nature of the transition or the observed patterns requires a more refined treatment. 
For instance, if the transition is discontinuous, the pattern might be very different from 
the linearly most unstable mode. This is why we perform a weakly nonlinear analysis 
close to the critical point. It is restricted to the diffusion limited regime (equations 12.111 
and l2.l2l2.lVt . 

In this analysis, nonlinear contributions to the evaporation rate are needed; they are 
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computed in appendix El 



J[h] 



1 



+Oih^). 



_a_ 

dx 

1 (ah\'^ ' 

2\dx/ 



(3.5) 



where Ti, is the Hilbert transform f appendix IB|| . Each successive correction to the base 
state is an integro-difFerential transform of the interface profile h; this comes from the 
non-local nature of the Laplace problem. 

We use a multi-scale expansion which is valid wh en the spatial Fourier spectrum of 
h{x,t) is concentrated around fee (see e.g lMannevill9.1990') . We look for an equation of 
evolution for the slowly varying function A{X, T) such that h{x, t) = A{X, T) exp{ikcx) + 
c.c. Formally, we use e as an expansion parameter. We assume that his a function of both 
the fast scales x, t and the slow scales X = ex, T = e^t. This choice for the slow scales 
is the natural one given that Q{k) is maximum at kc- We consider the neighbourhood of 
the marginal stability and we rescale the control parameter as 

nik,) = e^Ljik,). 

From the chain rule for differentiation, we make the replacements 

dx ~* dx + edx, dt ^ dt + e^dr- 

We also assume that h can be expanded as 

h{x,t) = 1 + ehi{x,t) + e^h2{x,t) H 

The procedure to obtain the amplitude equation (equation for A{X,T)) is quite stan- 
dard and is detailed in appendix |0 However the present case has the peculiarity that 
the expansion must be pursued up to order 6 as hi is found to vanish. This is due to 
the coupling between the evaporation rate and the fc = eigenmodes of the linearised 
evolution operator. 

Taking the rescalings X X/kc and T T/kc, neglecting terms of order k^ 
higher which is consistent with lubrication theory, the amplitude equation reads 



and 



dA 
df 



--a A 



3 d^A 



2 dX^ 
''-k^AlH 



2i 



dX^ 



d'^A 1 d^A 45 
2 ~ Y 

9 1 d 
dfr^dX 



A?A 



2i / Ji 



A 



dA 

'ax 



the control parameter being 



a = 2- 



Ca — Cuc 



Car 



3 Ma - Mac 
2 



Mar 



(3.6) 



(3.7) 



As far as we are aware of, this non-local kind of amplitude equation has not been derived 
before. The non-local terms are important for large amplitude patterns or for finite size 
systems. 

The solution A = of the amplitude equation becomes unstable when cr > (which 
is consistent with the linear stability). When > 0, there are stationary solutions of 
the form A{X, T) = Aq exp{iQX) with a = 45/8 \Ao\^ + 3/2 + 2Q^ + 1/2 Q*. They 
correspond to a stationary pattern with thickness fluctuation h{x,t) = ^oexp{i(fcc -|- 
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Q)x}, which is modulated around the critical wavenumber kc- Thus, as the prcfactor 
of is negative, the transition from the flat state {A = 0) to a state with height 

fluctuatio ns is a supercri t ical ( con tinuous) pitchfork bifurca tion in contrast with the 
studies of lVanHook et al\ l|l997|) or lThiele fc Knoblochl l)2004 on heated fluid layers. 



4. Discussion 

4.1. Transfer rate limited versus diffusion limited evaporation 

We obtained H2.9(l as a boundary condition prescribed at the interface combining both 
the transfer rate across the interface and the diffusion of the vapour in the gas phase. This 
boundary condition might be simplified in two distinct limits according to the values of 
the kinetic Peclet number and the thermal expansion number. These two limits have been 
used separately in the literat ure; however, it has be en overlooked that they fall within a 
general framework, although IMargerit et a,l\ l)200,ll) have treated vapour diffusion in the 
case of no evaporation. 

In the first limit, Pe\ <C min{x, k, fcx}, the evaporation process is limited by the trans- 
fer of molecules across the interface. The diffusion of the vapour can be ignored and the 
classical Hertz-Knudsen (|2.7|) relation gives the evaporation rate. It is worth noticing here 
that transfer-limited evaporation is always destabilising, if the free surface undergoes a 
small shape perturbation, liquid portions that are closer to (resp. farther from) the sub- 
strate evaporate faster (resp. slower) and the disturbance is amplified. The so-called one- 
sided models for e vaporating l avers jPr ospcrctti fc Plesset 1984; Burclbach et al. "19881; 
I Samid-Merzel et aL.1998.: Lvushnin et a/... 2002; .Margerit et aL 2003; Merkt fc BestehornI 
1200,31) correspond to this limit as they do not consider the dynamics of the gas phase and 
so discard vapour diffusion. 

In the second limit, Pey_ ^ max{x, A;}, the evaporation process is limited by the 
diffusion of the vapour in the gas phase. The gas phase is saturated in vapour immediately 
above the interface and the evaporation flux is given by the Fick relation 12.8|l . The 
electrostatic analogue of the Laplace problem for the vapour density, and the related 
sharp edge effect show that the evaporation rate is larger at bumps h ence diffusion 
limited ev aporation is stabilising. The study of evapora ting droplets by IPeegan et all 
|l999) and lCachile. Benichou. Poulard fc Cazaba'O l)2002|) fall within this approximation. 

It appears that this second limit is of large validity. Indeed, the ratio Pe-^/x which 
depends only on the nature of the fluid is the range 4-70 for the fluids considered here, 
whereas Pe\^ > 1 (a more restrictive condition) as soon as the thickness is larger than 
D/vth ~ 0.1 /im, which is the thickness scale determined by the diffusivity of vapour D 
and the thermal velocity t;th- 

4.2. Comparison with experiments on evaporating droplets 

We turn now to the experimental relevance of our analysis. Despite the number of the- 
oretical stud i es on evaporat ing thin fil ms, very few experiments hav e been conducted. 
iRedon et all l)l992|) . lKave"iK)Our et all l)200^ and lPoulard et all l)2003|) observed instabil- 
ities in the shape of evaporating droplets. In the first set of experiments ( Kavehpour et a/1 



|2QQ2; Redon et al. 1992), the spreading of drops of silicon oils was studied. As silicon oils 
have a very low volatility, evaporation rate is small and the s hear stress associate d with 
the spreading is important. In the second set of experiments dPoulard et a/Jl2003|) evap- 
o ration sets the v elocity scale. This is why we focus on this latter set. 

IPoulard et all l|2003j) studied receding evaporating drops of water and alkanes on 
smooth substrates. The data in tables 1 and 2 corresponds to this experimental situ- 
ation. For water, heptane and hexane the drop looses its axisymmetry and develops a 
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Ca 


Ma 




MadCa) 


Mac/Ma 


Water 


5.90 10"^ 


6.66 10- 


7 


8.59 lO"** 


12.9 


Nonane 


5.94 10"^ 


2.89 10" 


7 


8.60 lO"** 


30 


Octane 


1.4 10~** 


1.10 10" 


6 


1.52 10"^ 


13.8 


Heptane 


4.19 10"* 


3.76, 10" 


-6 


3.16 10-"^ 


8.4 


Hexane 


1.11 lO"'' 


1.33 10" 


5 


6.05 lO""^ 


4.55 



Table 3. Values of the control parameters (Ca, Ma) for different volatile fluids and 
comparison witfi the instabihty threshold in the diffusion hmited regime. 



regular wavy pattern near the contact line : the height fluctuates with a well-defined 
wavenumber. No instabilities are observed neither with octane nor nonane. We now com- 
pare with the stability analysis in the diffusion limited regime (13.411 . This comparison 
requires the choice of a thickness lengthscale; we retain the typical thickness of the un- 
stable region, h ^ 200 nm. Wc estimate the t ypical evaporation rate in this zone using 
J(r) ^ Jo/\/l - (Dccgan et ad ll999l) . R is the drop radius and r is the dis- 

tance to the drop axis. In the experiments, the radius R is of order 0.5 mm and the 
size of the unstable zone is i? — r ~ 5 /im. According to table 13 droplets of all flu- 
ids are stable. However, the geometry of a droplet is different from that of a constant 
thickness film analysed in this paper so we cannot conclude directly on the stability. We 
should only compare the relative values of the ratio Mad Ma between different fluids. 
This ratio is of order one meaning that the distance to the threshold is small. Moreover, 
the fluids can be sorted from the less stable to the most stable as: hexane, heptane, 
water, octane and nonane. This compares well with the experiments: festoon patterns 
have been observed with hexane, heptane and water and none with octane and nonane. 
However, caution should be t aken with water as e vaporating droplets of water have an 
anomalous retraction velocity ijPoulard et g^jEoO^ . Equation 1)3. 4|l gives as wavelengths 
Aneptane — 320 /zm and Awater — 600 jim which are an order of magnitude larger than 
the experimental Aneptane — 50 fim. and Awater — 30 jim. To summarize, the comparison 
with experiments is satisfying as we predict an instability threshold in contrast with the 
one-sided model (see section 3.2), and we flnd that the stability increases with the weight 
of the alkane. 

4.3. Main results 

In this paper, we constructed a two-sided model for evaporating thin liquid fllms. In 
order to predict quantitatively the evaporation rate, we have considered both transfer 
rate across the interf ace and diffusion of the vapour in the gas phase. The experiments of 
IPoulard et al\ ll2003|) motivated the study of the regime for which evaporation is diffusion- 
limited. In this context, the system describing the evolution of the height profile h couples 
the lubrication of the substrate by the thin film to the diffusion of its vapour in the 
gas phase. Using a linear stability analysis, we predicted an instability threshold and 
classified the stability for different liquids, in agreement with the experiments. The results 
of the linear stability analysis are not too far from Poulard et al. ( 2003) experiments on 
evaporating droplets. To push the comparison further, it would be interesting to perform 
experiments on extended flat films. 

Moreover, the diffusion equation confers a non-local aspect to the dynamics of the film 
which is found to persist in the amplitude equation H3.6(l establish ed within a weakly 
nonlin ear study. This property was unexpected and contrasts with iFriedrichs & Engell 
l|2003l) results on the Rosensweig instability of magnetic fiuids; it originates from the 
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presence of uniform profiles in the null-space of the linearised evolution operator (see 
appendix [C)l . 

This study suggests investigation of more complicated situations such as when both 
transfer rate and diffusion of the vapour are important (i.e. when Pek has a finite value), 
in the non-linear regime. As analytical extensions to the present study seem difficult, 
numerical computations of the complete two-sided model would allow the investigation 
of effects such as the non-stationary of the base state, the finite extension of the system, 
three dimensional patterns or the geometry of a droplet. Consequences of the non-local 
terms should also be interesting to further study. 

We are very grateful to Anne-Marie Cazabat and Christophe Poulard for getting us 
intersted in the stability of evaporating droplets. 



Appendix A. The evaporative flux. Electrostatic analogy. 

This appendix is devoted to the perturbative treatment of the Laplace problem asso- 
ciated with the diffusion limited regime 

V^p = (z > h{x)) 
p(x,h(x)) = 
§f(x,z = +w) = -l 

where h = h{x) is a given regular (bounded and differentiable) function. Precisely, the 
point is to compute — (n • ^)p\z=h{x)+ ■ 

A.l. Vortex sheet formalism 

Using electrostatic terminology, the problem is to find the electric field J immedi- 
ately above a deformed charged plane. Introducing a superficial charge distribution 
p = a{x)d{z — h(x)) {S{z — a) is a Dirac mass concentrated at the point z = a), one 
can write the integral representation: 



d£ = y I + {j^)^ dx being the arc-length element. 

At the interface, two boundary conditions are prescribed. To begin with, according to 
Gauss theorem, the field has a normal discontinuity 

df 

(J+ - J-) ■n = 2TTa— (A 2) 

dx 

where J'^ (resp. J~) stands for the field just above (resp. below) the interface z = 
h{x). Moreover, the tangential component has to vanish in order to fulfil the condition 
p{x, h{x)) = 0: 

J+t = J 1 = (A3) 

Now, setting z = h{x) in (|A l|l . one has to take the Cauchy principal value (denoted as 
PV) of the integral in order to get a well defined expression. This regularised integral is 
equal to the half-sum of and so that the conditions IjA 2IA 3p may be rewritten 
as: 
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PV 



X — x' 1 + x'Y 
with x') = and J+ ^ J+ ■ n ^ \\J+\\. 



-M{x) = 



(A 5) 



A. 2. Perturhative treatment 

The general relation (jA 5|) implicitly gives the superficial charge cr as a function of h. Since 
the inversion is not possible analytically, we assume that the deflection from the flat plane 
h is weak; we introduce a small parameter rj such as h is replaced by rjh (a possible choice 
is to set rj = sup^ \ h{x)\) and write a perturbative expansion cr = (t(°) +r/cr'^^^ +7y^cr(^) + . . .. 
The ?7 = contribution corresponds to the plane interface for which the field is uniform; 
the boundary condition J{x,z = +oo) = J^^-'e^ readily gives a (0) = j(")/7r. Solving 
l|X5|l up to 0(77'^), we find 



a[h] 



d/iV 
d^ / 



I— — 
dx^ dx 



n 



h-^n[h] 
dx 



(A 6) 



where Ti. is the Hilbert transform (see Appendix IB|) . 

The remaining step is simply to plug (|A 6|) into IjA 4p . Without having to compute 
further terms for a{x), we can express the evaporative fiux expansion up to C'(/i'*). 
Setting ry = 1, we find 



J[h] = 1 - ^n[h] 
dx 




d^^^d^^^d^^t^^] 




dx 



that is, the formula H3.5|) . 



Appendix B. Hilbert transform 

B.l. Definition and basic properties 

Given a bounded function f{x), we define the Hilbert transform with the usual conven- 
tions: 

n[f]ix) = - lim / dx'^^ 

TT e^O+y|^_^,|>^ X'-X 

where we have taken the Cauchy principal value (symmetric limit) at x' — x. Useful 
properties are commutation with linear differential operators and the inversion relation 
Ti.^^ = —H. With this definition, the Hilbert transform is not defined for constant func- 
tions. However, one can remove the divergence at infinity by taking the principal values 
both at X and at infinity; the result is then 7i[Cst] = 0. (Note that the inversion formula 
is not valid for constants.) 
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B.2. Hilhert transform and slow space varying amplitude 

In the weakly non- linear analysis, we have to compute quantities of the form H[A{ex) e*'^^], 
with < e <C 1. We want here to show that, up to a very good precision (for e sufficiently 
small), we have the relation 

7^[A(ex)e**^^] = A(ex)W[e*'=^] (fc ^ 0) (B 1) 

which means that the action of the Hilbert transform on a slowlv modulated Fourie r 
mode does not introduce, except if = 0, non-localities (see iFriedrichs fc Enge]ll20ni . 
In the next appendix, we make a substantial use of this property. 

Since A^{x) = A{ex) varies significantly only if x has a variation of order 1/e, the 
Fourier transform of A^ must be negligible outside of (— e, e). 

Let's first assume that the support of A^ (that is, the domain where it has non zero 
values) is included in (— e, e). Then, using 7i[e*'^^] — i sgn(fc) e**^^ (sgn(fc) — ±1 if ±/c > 0), 
we have 

n[A,{x)e''''']=ie'^'' j dfc'ie(fc')sgn(/c' + fc)e*'='^. 

Thus, IjB ip is proofed if we suppose |A:| > e, that is, if e is small enough. 

If = does not vanish outside of (— e,e), one can show that, under the same 
hypothesis on k, to the first order in e, the error on ETl) is of the order of /j7, |i(A;)|dfc. 



Appendix C. Weakly non linear analysis 

Here we detail the weakly non linear analysis leading to the amplitude equation (|3.6(l . 
We consider the neighbourhood of the stability limit, so that we rescale the control 
parameter according to 

n{K)^e^u{K). (CI) 

Plugging (|3.5() into the lubrication equation l|2.11() , we obtain a closed integro-differential 
equation for the height profile h. We assume that ft, is a function of x,X,t,T (fast and 
slow scales) and admits the expansion h = e/i^^^ + e^h^"^^ + ■ • • . The derivatives are 
substituted according to 

d^^d^ + edx, dt^dt + e^dr. (C 2) 

It is convenient to separate the evolution operator into its linear contributions C and 
non-linear ones A/". In particular, applying transformations (|C 2p leads to the expansion 

^ - iCa'^''- 

Note that the null-space of Cc contains slow space varying height profiles (i.e. functions 
of X). We now proceed to the solution order by order. 
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Order e^; 

We have simply: 

L^h^^^^Q. (C3) 

Using l|B the solution is 

= (An(X,r)e*'=^^+ c.c)+Aio(X,T) 
kc being the critical wavenumber given by the linear stability analysis. 

Order e^; 

The equation has the form 

C,h^^^ ^-C^^)h^^^ -N^-'^h^^^) (C4) 

The non-linear term contains a fc = mode. As the right-hand side of (jC 4|) must be 
orthogonal to the null space of Cc, it imphes /i'^^ — 0. Thus, Cch^"^^ = 0, hence h'^^^ = 
{A21 (X, T) e'^^c- + c.c) + A20 (X, T) . 

Order : 

AA(3) = as = 0: 

C,h('^ = -C^'^h^'\ (C5) 

Again, right-hand side of IjC 5|l has to be orthogonal to the null space of Cc, so we have 
A20 = 0. Hence, /i^^) = (A3i(X, T) e''^'^^ + cc) + A3o{X,T). 
Order : 

The non linear term is: 

A/-W(/i(2)) ^ [f^Al^e''"'^'' + cc) - klA2{Mi 
with {A21 is the complex conjugate of A21) 

So, introducing a = {~-^kl + ^(1 + S/cc)), we have from lfU6)l 

-axW[A30] =fc2|A21P, (C7) 

-w(fcc)A2i + ad\A2i + 9tA2i = 0. (C 8) 

The solution to IjC 6|l is 

hi = ( A42 [X, T) e2^'==- + cc) + {Ail {X, T) e^'=^" + cc) + A40 T) 
with Ai2 = I^Ali. 

Order : 

LM^^ = -£(3)/,(2) _£(2)^(3) _^(l);^(4) __^(5)(;^(2)^;^(3))^ 9) 

It is convenient to decompose the non-linear term: 

AA(5)(/j(2)) ^ ^^^^ ^2^K. ^ ^ ^^^^ ^ Cc) + A/'so, 



where 



and with 
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AA52 = /5^3i^2i + igsdxAl^ 

= i5^30^21 

A4o = -^2^21^;^ - ihA2idxA^ - kcH [A^dxA2i] + c.c 

/5(fcc) = ^ (-8fc,^ - 2k'^ - Akl) + ^ + kl 
9M = ^{2k^ + 2kl + &kl)-^^-^-l 

Imposing again that the restriction of the right-hand side of IjC 9|) to the nuU space of Cc 
has to vanish and using (|C 7p . we get 

il3d\A2l + - C^(fcc) + ^t} v431 + J5^21 / W[|A2l|2] = (C 10) 

and 

dxn[Aio] = (-C^(fce) + 5t)^30 + ^^i^3o +A4o. (Cll) 

The solution at this order is 

h5 = (A52 (X, T) e2^'==^ + c.c) + ( Ai (X, T) e*'=^^ + c.c) + Ao (X, T) 

with A52 = (^^52 + if^'(2fcc)9x^42). 

Even if (ICIfl is non-hnear in A21, it does not give the nature of the bifurcation. This 
is why we carry on computations to next order. 

Order : 

CM""^ = -£(4);^(2) „^(3)^(3) _^(2)/^(4) _^(l)/^(5) _A/-(6)(/j(2)^/j(3)^^(4))^ 12) 

We only need the part of TV-^-* of wavenumber k^. 

J^&l = /e ^2 1^42 + 56 ^21 9x^^30 + ije ^21 9x^30 

+4(^21^40 + A31A30) + mG\A2ifA2i + ineA3odxA2i 

with 



feikc) 


^ Ma I 1,2 


9 1,3 


/ ^ Ca 




gdkc) 


_ Ma (1 1,3 1 


ll,2^ 


1 ~h 




jdkc) 


= 2M^{k! + , 


kl + 


'^^) Ca 




(eikc) 


_ Ma (3 1.3 
~ Ca \2'^c 


+ k'. 


/ ^ Ca 


- jhikc) 


nieikc) 


_ Ma (7 1.5 
~ Ca U'^c 


1 l2 
2'^c 


3l4 


¥1) 


ne{kc) 


_ Ma (9 1.2 I 


2k,) 


ikl 
Ca 





7k: 



After solving equations H(' 7I(' 11(1 for A30 and ^40, we use the solvability condition that 
the right-hand side of ((C 12(1 is orthogonal to the null space of Cc and we obtain an 
equation for A41. Introducing 



A = eAii + e^A2i + e^Azi + e^A^i 
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the last equation for An can be re-summed with equations (|C 8IC 10(1 for A21 and A31. 
Using the inverse transformations of (|C IIC 2\ . we finally obtain the amplitude equation 
for A{X,T): 



i~nik,) + adj. + dT)A + i(3d^^A + —d%A + ^l^pA 
- WU\'']A + in^klidxA) I n[\A\^] 
+£e!^k'',A J j (-O(fcc) + dT + dxn)\A\'' + 2ik,A J H [AdxA\ | = 0. 



with 



n{2kc 



:/6/2 



J Ma 
2C^ 



2kr I ii 



Taking the limit of small fcc with the help of the rescaling dx k^dx and Ot — > kcdx 
leads to 



d_ 

dT 



3 9^ \ 
2 9X2 j 



, „ 1 45, ,,2 , ^, .,21 . 



—.K{A \—-a + —H\\A\^ + 2^A H 



dT 



dX 



A 



dA 
dX 



(C13) 



hence the simplification ((3.6|l . Note that l(C 13|l is valid only if A{X,T) vanishes at X = 
±00, so that Hilbert transforms Ti. are well defined. In contrast with standard weakly 
non-linear analysis, this equation is non-local. 



REFERENCES 

Benard, H. 1900 Les tourbillons cellulaires dans une nappe liquide. Rev. Gen. Sci. Pures Appl. 
11, 1261-1276. 

Berg, J. C, Boudart, M. & Acrivos, A. 1966 Natural convection in pools of evaporating 
liquids. J. Fluid Mech. 24, 721-735. 

BuRELBACH, J. P., Bankoff, S. G. & Davis, S. H. 1988 Nonlinear stability of evaporat- 
ing/condensing liquid film. J. Fluid Mech. 195, 463-494. 

Cachile, M., Benichou, O., Poulard, C. & Cazabat, A. M. 2002 Evaporating droplets. 
Langmuir 18 (21), 8070-8078. 

Davis, S. H. 1987 Thermocapillary instabilities. Annu. Rev. Fluid Mech. 19, 403-435. 

Deegan, R. D., Bakajin, O., Dupont, T. F., Huber, G., Nagel, S. R. & Witten, T. A. 
1997 Capillary flow as the cause of ring stains from dried liquid drops. Nature 389, 827-829. 

Deegan, R. D., Bakajin, O., Dupont, T. F., Huber, G., Nagel, S. R. & Witten, T. A. 
1999 Contact line deposit in an evaporating drop. Rhys. Rev. E 62 (1), 217-239. 

Elbaum, M. & LiPSON, S. G. 1994 How Does a Thin Wetted Film Dry Up? Rhys. Rev. Lett. 
72 (22). 

Fanton, X. & Cazabat, A. M. 1998 Spreading and instabilities induced by a soluted 

Marangoni effect. Langmuir 14, 2554-2561. 
Fournier, J. B. & Cazabat, A. M. 1992 Tears of wine. Europhys. Lett. 20, 517-524. 
Friedrichs, R. & Engel, A. 2003 Nonlinear analysis of the Rosensweig instability. Europhys. 

Lett. 63 (6), 826-832. 

GOLOVIN, A. A., Nepomnyashchy, a. A. & Pismen, L. M. 1997 Nonlinear evolution and 
secondary instabilities of Marangoni convection in a liquid-gas system with deformable 
interface. J. Fluid Mech. 341, 317-341. 

Hegseth, ,J. J., Rashidnia, N. & Chai, A. 1996 Natural convection in droplet evaporation. 
Rhys. Rev. E 54 (2). 



Evaporation instability in a thin film 



19 



Hosoi, A. E. & Bush, W. M. 2001 Evaporative instabilities in climbing films. J. Fluid Mech. 
442, 217-239. 

Kalliadasis, S., Kiyashko, a. & Demekhin, E. A. 2003 Marangoni instability of a thin liquid 
film heated from below by a local heat source. J. Fluid Mech. 475, 377-408. 

Kavehpour, P., OvRYN, B. & McKlNLEY, G. H. 2002 Evaporativcly-drivcn Marangoni insta- 
bilities of volatile liquid spreading on thermally conductive substrates. Colloid Surface A. 
206, 409-423. 

Langer, J. S. 1980 Pattern formation in crystal growth. Rev. Mod. Phys. 52 (1), 1-28. 

Lide, S. R. 1995 Handbook of Chemistry and Physics, 75th edn. CRC Press. 

Lyushnin, a. v., Golovin, A. A. & Pismen, L. M. 2002 Fingering instability of thin evapo- 
rating liquid films. Phys. Rev. E 65. 

Mancini, H. & Maza, D. 2003 Pattern formation without heating in an evaporative convection 
experiment. Europhys. Lett. 66 (6), 1-7. 

Manneville, p. 1990 Dissipative structures and weak turbulence. Academic Press. 

Marangoni, C. 1865 On the expansion of a drop of liquid floating on the surface of another 
liqmd. Pavia: Tipographia dei fratelli Fusi. 

Margerit, J., COLINET, P., Lebon, G., Iorio, C. S. & Legros, J. C. 2003 Interfacial 
nonequilibrium and Benard-Marangoni instability of a liquid vapor system. Phys. Rev. E 
68. 

Marra, J. & Huethorst, J. A. M. 1991 Physical Principles of Marangoni Drying. Langmuir 
7, 2748-2755. 

Matar, O. K. & Craster, R. V. 2001 Models for Marangoni Drying. Phys. Fluids 13, 1869- 
1883. 

Merkt, D. & Besteiiorn, M. 2003 Benard-Marangoni convection in a strongly evaporating 

fluid. Physica D 185, 196-208. 
MiLADiNOVA, S., Slavtchev, S., Lebon, G. & Legros, J. C. 2002 Long-wave instabilities of 

non-uniformly heated falling films. J. Fluid Mech. 453, 153-175. 
Nguyen, V. X. & Stebe, K. J. 2002 Patterning of small particles by a surfactant-enhanced 

Marangoni-Benard instability. Phys. Rev. E 88 (16). 
O'Brien, S. B. G. M. 1993 On Marangoni drying: nonlinear kinematic waves in a thin film. J. 

Fluid Mech. 254, 649-670. 
Oron, a., Davis, S. H. & Bankoff, S. G. 1997 Long-scale evolution of thin hquid films. Rev. 

Mod. Phys. 69 (3), 931-980. 
Poulard, C., Benichou, O. & Cazabat, A. M. 2003 Freely receding evaporating droplets. 

Langmuir 19 (21), 8828-8834. 
Prosperetti, a. & Plesset, M. S. 1984 The stability of an evaporating liquid surface. Phys. 

Fluid 27, 1590-1602. 

Redon, C., Brochard-Wyart, F. & Rondelez, F. 1992 Festoon instabilities of shghtly 

volatile liquids during spreading. J. Phys. II 2 (9), 1671-1676. 
Samid-Merzel, N., LiPSON, S. G. & Tanniiauser, D. S. 1998 Pattern formation in drying 

water films. Phys. Rev. E 57 (3). 
SCHATZ, M. F. & Neitzel, G. p. 2001 Experiments on thermocapillary instabilities. Annu. 

Rev. Fluid Mech. 33, 93-129. 
Thiele, U. & Knobloch, E. 2004 Thin liquid films on a slightly inchned heated plate. Physica 

D 190, 213-248. 

Thomson, .J. 1855 On certain curious motions observable at the surfaces of wine and other 
alcoholic liquors. Phil. Mag. 10, 330. 

VanHook, S., Schatz, M., Schatz, M. F., Swift, J. B., McCormick, W. D. & Swinney, 
H. L. 1997 Long-wavelength surface-length-driven Benard convection: experiment and the- 
ory. J. Fluid Mech. 345, 45-78. 

Vuilleumier, R., Ego, V., Neltner, L. & Cazabat, A. M. 1995 Tears of wine: the stationary 
state. Langmuir 11, 4117-4121. 

Yeo, L. Y., Craster, R. V. & Matar, O. K. 2003 Marangoni instability of a thin liquid film 
on a locally heated horizontal wall. Phys. Rev. E 67. 



